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Summary 

Two perturbation techniques are applied to two singular perturbation prob- 
lems in heat transfer to obtain uniformly valid solutions which can serve as 
benchmarks for the numerical techniques: finite-difference and finite-element 
techniques. In the first problem, the method of strained parameters coupled 
with the application of a solvability condition is used to obtain a uniform 
solution for the problem of unsteady heat conduction in a long nearly circular 
cylinder. In the second problem, the method of matched as^ptotic expansion 
coupled with Van Dyke's matching principle is used to obtain a uniform solution 
for the problem of one dimensional conduction-convection heat transfer of a 
uniform fluid flow. 


I. Introduction 


The main purpose of this paper is to demonstrate the capabilities of the 
perturbation techniques in developing approximate closed- form solutions for 
heat transfer problems involving difficulties which preclude their solutions 
exactly or require resorting to computational techniques such as finite- 
difference, finite-element, and panel techniques. These difficulties may be 
due to nonlinear governing equations, equations with variable coefficients, 
nonlinear boundaries, and existance of boundary layers near portions of the 
boundaries . 

Although computational techniques in the areas of fluid dynamics, heat 
transfer, and structures are rapidly advancing and are capable of developing 
excellent solutions for realistic problems, one always needs bench-mark 
solutions, if experimental data are not available, to check the developed com- 
puter code or to check the accuracy of the computed results. In this regard, 
a closed form perturbation solution for a simplified problem which retains the 
same difficulties (weakly nonlinear equations and boundary conditions and 
weakly irregular boundaries) can best serve this purpose. 

Among the perturbation techniques, the strai ghtforward expansion in terms 
of a parameter in the problem leads to satisfactory results if one is dealing 
with a regular perturbation problem or if its region of nonuniformity is 
avoided (ref. 1 and 2). However, for singular perturbation problems, the 
straightforward expansions yield nonuniform solutions and one has to use 
other perturbation techniques to obtain uniform solutions (ref. 3-6). Infinite 
domains in a problem, a small parameter multiplying the highest derivative of the 
governing equation, type change of a partial differential equation, and 
existence of singularities are some of the sources of nonuniformities of 
strai ghtforv/ard solutions. 
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In this paper, two applications in the area of heat transfer are con- 
sidered and closed-form uniform perturbation solutions are developed. In the 
first application, the problem of unsteady heat conduction in a long 
nearly circular cylinder is considered and a straightforward solution is 
shown to breakdown. The method of strained parameters coupled with the 
application of a solvability condition is used to develop a uniform 
solution. 

Similar problems in the areas of duct acoustics (ref. 7) and vibrations 
(ref. 3) were considered where the methods of multiple scales and the method 
of strained parameters were used, respectively. 

In the second application, the problem of one dimensional conduction-con- 
vection heat transfer of a uniform fluid flow in a single channel is considered. 
For small ratios of conduction to convection heat transfer, the problem 
is shown to possess a thermal boundary layer where large temperature gra- 
dients exist. The method of matched asymptotic expansion coupled with Van 
Dyke's matching principle is used to develop a uniform solution. This problem 
was considered in reference 8 for single and merging flows by using the finite- 
element technique. Steady tv/o-dimensional problems with different locations 
of the boundary layer can be found in reference 6. 
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II. Un steady Heat Conduction in a Long Nearly Circular Cylinder 


We consider the two-dimensional unsteady heat conduction in a long cylin- 
der whose cross sectional area is nearly circular. Initially, the cylinder is 
at temperature g(r*,(f>) and at any later time the surface is kept at zero tem- 
perature. The radius of the cylinder is expressed as 

2tt 

r* = R + a f((|)) where / f((|)) d <|) * 0 and a « R (1) 


Dimensionless quantities are introduced by using the mean radius of the cylin 
der R, the characteristic temperature T , and the time R^/a (a is the thermal 

diffusivity) as reference quantities. The dimensionless form of the problem 
is given by 


1 ~ 1 - 

0 +-0+-e =0 

rr r r ^2 ^cfxj) 


0(rQ.4>,t) = 0 


on r. 


( 2 ) 


2tt 

+ e f((j)) and / f((j)) d<{) = 0 (3) 

0 


e(r,(j), o) = g(r,(f)) 


(4) 


The parameter e is a small quantity characterizing the small deviation of the 
cross sectional area from the circular shape. The temporal Variation is sepa- 
rated by assuming a solution of the form 

0(r,(j),t) = e(r,<p) (5) 

Substituting equation (5) into equations (2) and (3), we obtain 


rr 


1 

r 


H 


+ B e = 0 


( 6 ) 


=0 on r^ = 1 + e f((j>) (7) 

Equation (6) is the Helmholtz equation. Although equation (6) is linear, the 
problem is not separable because the boundary condition, in the present form, 
is not separable. Since e is a small parameter, one can expand eCr,4);e) in 
the form of a power series in terms of c as follows 

0(r,4);e) = 0^ (r,4.) + e 0^ (r,4>) + (8) 
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In equation (8), only two terms are considered and hence a first-order solu- 
tion is intended. Since e appears in the argument ofe, equation (7), one 
needs to extract e from this argument so that the process of equating coeffi- 
cients of like powers of e can correctly be accomplished. Therefore, the 
boundary condition of equation (11) is expanded around r = 1 using a Taylor- 
series expansion. This process is well known as the "transfer of the boun- 
dary condition." Thus, we get 


0(rQ,<|)) = e(l,<|)) + e f((|)) 0^(1, (|)) + = 0 (9) 

Substituting equation (8) into equations (6) and (9) and equating coefficients 
of like powers of e, we obtain the following two sets of problems: 

0(e°) - Problem 


0 + “ 9 + ~n 0 ,, + 3^0 = 0 

orr r or ^2 o4>(|) o 


0q(1,4>) = 0 


0(g) - Problem 


®lrr ^ r ^Ir ^ ^2 ^l(f>(t> + 3^9i = 0 
0l(l,<f)) = - f((|)) 0Q^ (l,<f>) 


( 10 ) 

(11) 


( 12 ) 

(13) 


In the perturbation expansion used above, we note that only the dependent 
variable 0 is expanded in terms of the small parameter e. Such a perturba- 
tion method is called a "straightforward-perturbation method." Straightfor- 
ward expansions break down when we deal with singular perturbation problems. 
Next, it is shown that the straightforward expansion breaks down for this 
problem. 

The solution of the 0(e°) problem is obtained by using the method of 
separation of variables. The solution is found as 


o n nm nm nm 


(14) 


where k^ (=3) are the zeros of the J (3) = 0, and A are complex constants, 
nm n nm 

Setting 3^ = k^^^ in equation (12) and subtituting equation (14) into 
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equation (13), we get 


11 ? 

9, + - 0, + -o 0 + k 0, = 0 

>■ 1# ™ ’ 

e.d,*) = - k J'(k ) t{J>) (A e*"* +S e’''’*) 

1 nm n' nm' ' nm nm 


(15) 

(16) 


The solution of the 0(e) problem is obtained by expanding each of 0, and 
f(<j>) in a Fourier series as 


0l(r,t()) =_? G^(r) e 


it(j) 


(17) 

(18) 


1 


2tt 


where fp = ^ ■/’ f('f>) and f q = 0 according to the condition in equa- 
tion (3). Substituting equation (17) and (18) into equations (15) and (16), 
multiplying the results by expC-iscj)) and integrating from 0 to 2 tt, we get 


Q"+ 1 G' + (k - -^-) G = 0 

s r s ' nm 2' s 

r 

65(1) *^n ^*^nm^ ^^nm ^s-n ^nm ^s+n^ 


(19) 

(20) 


If s ji* n, equations (19) and (20) have a unique solution since the only solu- 
tion for the corresponding homogeneous problem is the trivial solution. If 
s = n , we obtain 


K' + l K + Ck„„ 

n r n nm 


G„ = 0 


( 21 ) 


■ ■ ^n, J;; 


( 22 ) 


we note that the corresponding homogeneous problem has a nontrivial solution 
and hence the inhomogeneous problem will not have a solution unless a solva- 
bility condition is satisfied. However, a solvability condition does not 
exist and hence there is no solution to equations (21) and (22). Therefore, 
the solution of the 0(e) problem breaks down. 
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The reason behind this trouble is due to the straightforward method used 
here. To obtain a uniform solution, we use the method of strained para- 
meters. In this method, we expand 3, in addition to the expansion given 
by equation (8), as follows 


3 = 3q + e 3^ + — 


( 23 ) 


where 3^ is to be determined in the course of the solution. Substituting 

equations (8) and (23) into equation (6), substituting equation (8) into 
equation (9), and equating coefficients of like powers of e, we again obtain two 
sets of problems. The 0(e®) problem is the same as that given by equations 
(10) and (11) with the exception of replacing 3^ by 3^. The 0(e) problem 
is given by 


Irr 


+ - e, 

r 1 r 




+ b; 


6i = 


- 2 B7i 


(24) 


and by equation (13). 

Substituting equation (14) into equation (24) and replacing 3^ by 
we get 


0, + — 6, + — 9 Oil. kf 0-| - 

Irr r Ir ^2 Icjxf) nm 1 


-2k 3, J (k r)(A + A e''''^'^’) 

nm 1 n nm nm nm ' 


(25) 


The boundary condition is still given by equation (16). 

The solution of equations (25) and (16) is obtained by substituting 
equations (17) and (18) into equations (25) and (16), multiplying the result 

by exp(-is(f>) and integrating from 0 to 2 tt. Again, we obtain two cases 

corresponding to s jf^ n and s = n. The former case is the same as that given 

by equations (19) and (20) in which a unique solution exists. In the latter 

case, we have 




= - 2 k„^ 3l Kn. 

2 n nm I nm n nm 


(26) 


' - ''nm \m ''n'f'nm) 


( 22 ) 


As we mentioned before, the problem given by equations (26) and (22) has a 
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solution if and only if a solvability condition exists. To obtain the solvability 
condition, we writd equation (26) in the self-adjoint form by multiplying both 
sides of the equation by r. Next, we multiply both sides by the adjoint u(r) 
and integrate the result by parts over the range of r. Thus, we get 


1 

/ G 
0 


n 


[(ru-)- + (k 


2 

nm 


r 


n2 1 

5- )u] dr + [r u G;; - r u" Gj 


^ *^nm ^1 ^nm 


1 


(27) 


The adjoint equation is obtained by setting the coefficient of G^^ to zero. 

Thus, we obtain 

2 

(r u")^ + ■ r ^ u = 0 (28) 

The adjoint boundary conditions are obtained by choosing 

u(l) = 0 and u(0) < <» (29) 


in equation (27), Equation (27) reduces to 


(30) 


The solution of the adjoint problem, equations (28) and (29), is given by 

u(r) = J„(k ^r) 
n nm 

Substituting equations (22) and (31) into equation (30), and performing the 
integration on the right hand side, we obtain the equation defining 3-| as 


I nm nm nm 


(32) 


since fo„, and are complex constants, we assume 
2n nm nm 

f - h J'^2n 
^2n ■ ^2n ® 


A 1 ^ T ^ a o ^ 

^nm " 2 ®nm ® ’ ^nm " 2 ®nm ® 


(33) 
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Substituting equation (33) into equation (32), and equating the real and imaginary 
parts, we get 


^nm "" 2 ^2n 


or 



7t) 


(34) 


~*^nm *^2n 


or 


*^nm ^2n 


(35) 


Substituting equation (34) into equation (33) and substituting the result into 
equation (4), we get 


e 


( 1 ) 




(36) 


and 

' ^nm <" (37) 

Substituting equation (35) into equation (23), we get 


- e k b„ 
nm nm 2n 


6^^^ = k + e k b„ 
nm nm 2n 


Substituting equation (36) into equation (8), substituting this result and 
equation (38) into equation (5), repeating the same process with equations 
(37) and (39), and forming a linear combination of the two, we obtain 


(38) 

(39) 


9 (r.*.t) • 

+ sinh ♦ + Y2„/2) Exp[- k^(k^ + 2e + — - (40) 

The constants a ^ and c are found from the initial condition, 
nm nm 
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One Dimensiona l Conduction-Convec ti on fleat-Transfer in a Uniform Flow 

We consider the one dimensional convective diffusion equation modified by 
a convective surface loss term. This equation has been used to model the far- 
field behavior of thermal regime for single and merging fluid flows. Here, 
only the problem of single channel flows with specified temperatures at the 
upstream and downstream boundaries is considered. The governing equation of 
the average temperature T(x*) is given by 


-k A T"" + p C A u T' + h p(T - Tg) = 0 (1) 

The boundary conditions are 

T(0) = T^ (2) 

T(L) = (3) 

In equations (1) and (2), k is the coefficient of thermal conductivity in the 
flow direction, A is the flow cross-sectional area, p is the fluid density, C 
is the fluid specific heat, u is the flow average velocity, h is the convec- 
tion heat exchange coefficient, p is the convection perimeter, T^ is the con- 
vection exchange temperature, and L is the channel length. 

Dimensionless parameters are introduced by using the pipe length L, and the 
temperature difference T, -T as reference quantities. The dimensionless form 
of the problem is given oy 

-n/Pg) r' + y" + ^ ° 

y(0) = 1 , y(l) = 02 (5) 

where = pCuL/k is the Peclet number, = hd/k is the Nusselt number, 

d = 4A/b is the hydraulic diameter, and y(x) = [T(x) - Tg]/(T.| -T^) . For 

large Peclet numbers (small ratio of conduction to convection heat transfer) 
we let 1/Pg = e, where e is now a small parameter. Moreover, we assume 

N^/P^ = b = 0(1). 

Thus, we obtain the following problem describing the spatial variation 
of the temperature: 


e - y" - b y = 0 

(6) 

y(o) = 1 , y(l) = 02 

(7) 
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Equations (6) and (7) describe a typical boundary-layer problem where a 
small parameter multiplies the highest derivative. This problem can success- 
fully be treated by using any of the several suitable perturbation techniques, 
namely, the method of matched asymptotic expansion, the method of multiple 
scales, and the method of composite expansions, among others. In this paper, 
we develop a uniformly valid solution by using the method of matched asymptotic 
expansion . 

Since the coefficient of y' is negative in the interval 0 £ x j< 1 , the 
boundary layer exists at the boundary x = l. Next, we develop outer and inner 
solutions and match them by using Van Dyke's matching principle to obtain a 
composite solution which is uniformly valid everywhere. 

Outer Solution : 

The outer solution y°(x;e) is expressed in the form 


y°(x;e) = y^(x) + e y^(x) + (8) 

Dropping the second boundary condition of equation (7), substituting equation 
(8) into equation (6) and into the first boundary condition of equation (7), 
and equating coefficients of like powers of e, we obtain the following prob- 


1 ems : 

0(£°) Problem 

+ b = 0 (9) 

y„(o) = 1 (10) 

O(e^) Problem 

y{ + b y^ = y^" (11 ) 

y^(0) = 0 (12) 

The outer solution of equations (9)-(12) is given by 

y° = e~^^(l + e b^ x) + O(e^) (13) 
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Inner Solution: 


To develop an inner solution valid near the boundary layer, we drop the 
boundary condition at x = 0 and introduce the following stretching transforma- 
tion for the inner variable 

5 = ’r 

(14) 

In terms of the inner variable, the inner 

solution y^(5;e) is governed by 

^ i =0 

(15) 

y'(o) = 02 

(16) 

Next, we expand the inner solution in the 

form 

y^5;e) = Yq(5) + E Y^(C) + 

(17) 

Substituting equation (17) into equations 
cients of like powers of e, we obtain the 

(15) and (16) and equating coeffi- 
fol lowing problems: 

0(e°) Problem 


Y" = 0 
0 0 

(18) 

Vo(0) = 02 

(19) 

O(e^) Problem 


'i'' + y: = b Y 
1 1 0 

(20) 

O 

II 

o 

>- 

(21) 


The inner solution of equations (18)-(21) is given by 

y’ = Go + A^(l - e"^) + e{b C[02 + A^Cl + e'^)] + A^(l - e"^)} + O(e^) 

( 22 ) 

We note that the outer solution is completely known while the inner solu- 
tion contains the two unknown coefficients A^ and A^ . They are determined by 
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applying Van Dyke's Hatching principle to the inner and outer solutions. The 
principle states that: 

the m-term inner expansion of (the n-term outer expansion) = the n-term outer 
expansion of (the m-term inner expansion) 

where m and n are any two integers. To find the left hand side, we write the 
n-term outer expansion in terms of the inner variable, equation (14), expand 
the functions of e keeping the inner variable fixed, and keep m-terms of the 
resulting expansions. An opposite procedure is applied to the right hand 
side. For m = n = 2, we obtain 

Aq = e"*^ - 02. = b^ e"^ (23) 

Substituting equation (23) into equation (22), we obtain the inner solution as 

= 02 + (e“^ - 92)(1 - e"^) + e{b ^[02 + (e‘*^ - 02)C' + e“^)] + b^e~^(l-e'^)} 

+ O(e^) (24) 


Next, we express the composite expansion in the form 

y"^ = y° +y' - (y°)' (25) 


Substituting equations (13) and (24) and the result of the left hand side of 
Van-Dyke's matching principle in equation (25), we obtain the composite solu- 
tion as 

- ( 1 ) /e 

y^ = (1 + e b^x)e"^^ + [(l-b+bx)( 02 -e"‘^) - e bV^^® + O(e^) • (26) 

It should be noted that an exact solution exists for equations (6) and (7) 


which is 

given by 



y = 

iDq X 

[(e ^ - 02 )e ^ 

itIqX itIq 

- (e ^ - 02 )e 2 ]/(e ^ - e ^) 

(27) 

where m^ 

^2 = n ± e(l = - 

1 e b)^^^]/2 e 

(28) 

The 

solution of the 

unsteady flow problem described by 


-k A T * 

X*) 

,* + P C A U T^* 

+ h p(T-Tg) + p C A T^* = 0 

(29) 
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(30) 


T(0,t*) = . T(L,t*) = T2 , T(x*,t*) = 

is an easy extension to the steady solution given by equation (26). The prob- 
lem is divided into steady and transient problems. The steady problem with 
the inhomogeneous boundary conditions is already considered. Upon assuming an 
exponential time-decay solution, the transient problem with homogeneous boun- 
dary conditions reduces to an eigenvalue problem coupled with a boundary 
layer. A uniform solution of the problem is obtained by using the method of 
matched asymptotic expansions (as shown before) coupled with strained eigen- 
values which are expressed as a power series of e. 


IV. Concluding' Remarks 

The methods of strained parameters and matched asymptotic expansions are 
successfully used to obtain closed-form perturbation solutions for heat- 
transfer problems with irregular boundaries and with boundary layers, respec- 
tively. The techniques given here are applicable to a large class of similar 
problems where various geometrical shapes and three-dimensional dependence are 
considered. Nevertheless, such solutions serve as bench marks for the compu- 
tational techniques. Moreover, useful solutions can be obtained by combining 
a computational technique with a perturbation technique in a certain problem. 
The computational technique is used to solve the zeroth-order problem while 
the perturbation technique is used to solve the higher-order problems. 
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